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The level densities and y-ray strength functions of 205-208p^ hssa been measured with the Oslo method, 
utilizing the (^He, ■'He'y) and (■'He, ay) reactions on the target nuclei ^^''Pb and ^^**Pb. The extracted level 
densities are consistent with known discrete levels at low excitation energies. The entropies and temperatures 
in the micro-canonical ensemble have been deduced from the experimental level density. An average entropy 
difference of A5 ~ 0.9A;g has been observed between ^"^Pb and ^'^^Pb. The y-ray strength functions in 205-208 
have been extracted and compared with two models; however, none of them describe the data adequately. Inter- 
mediate structures have been observed in the y-ray strength functions for y-ray energies below neutron threshold 
in all the analyzed Pb nuclei. These structures become less pronounced while moving from the doubly-magic 
nucleus ^O^Pb to ^'Wpb. 

PACS numbers: 21.10.Ma, 21.10.Pc. 25.55.Hp, 27.80.+W 



I. INTRODUCTION 

Tracing average nuclear properties such as nuclear level 
density and y-ray strength function in the quasi-continuum re- 
gion are of particular importance. Both level density and y-ray 
strength function are inputs in statistical model calculations of 
compound nuclear reactions and the subsequent decay of the 
compound system. These calculations are important for many 
aspects of nuclear structure studies, e.g., fission hindrance in 
heavy nuclei, giant resonances built on excited states, yields 
of evaporation residues to populate certain exotic nuclei and 
production of heavy elements in stellar processes. The level 
density is also an essential quantity for determining thermody- 
namic properties of nuclei, such as entropy and temperature - 
quantities that describe the many-particle behavior of the sys- 
tem. 

The level density is defined as the number of levels per unit 
of excitation energy. It can be obtained experimentally from 
different methods such as counting of resonances following 
neutron capture 1 1 ] and modeling of particle evaporation spec- 
tra from compound nucleus reactions The y-ray strength 
function characterizes the nuclear electromagnetic response. 
The concept of y-ray strength function was introduced in the 
work of Blatt and Weisskopf [3]. Most of the experimental 
information on the y-ray strength function has been obtained 
from the study of photonuclear cross-sections for high energy 
y transitions {Ey ~ 10 — 20 MeV) |4j]; however, at low y ener- 
gies experimental data are scarce. 

The Oslo cyclotron group has developed a method ^ to 
isolate the first y-rays emitted in all decay cascades at var- 
ious initial excitation energies. The energy distribution of 
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these primary (or first-generation) y-rays provides informa- 
tion on the level density and the y-ray strength function. The 
Oslo method has previously been applied for rare earth nu- 
clei, obtaining information on the level density and average 
y-decay properties in this mass region [6, 7]. Recently, the 
method has been extended to other mass regions like Fe, Mo, 
V and Sc S H [S IBi JJJ. In order to check the validity of 
the method in cases where nuclei have low level density and 
strong Porter- Thomas fluctuations [13], the ^^^^^Si nuclei were 
studied with gratifying results [14]. These achievements have 
encouraged us to apply the method on closed shell 205-208p^ 
nuclei, where the decay properties are less statistical due to 
shell effects and more dominated by single-particle selection 
rules. 

Based on the shell model, the nucleon configuration of 
lead isotopes is such that the protons fill the /111/2 shell, and 
the valence neutrons reside in the /13/2 shell. The chain of 
A — 205 — 208 lead isotopes in general, and the doubly-magic 
208p^ nucleus with Z = 82 and — 126 in particular, are in- 
teresting nuclei for studying the nuclear structure and y-decay 
properties at and in the vicinity of a two major shell closures. 

In the following sections a brief outline of the experimental 
method and data analysis are given. The level densities, ther- 
modynamic properties and the y-ray strength functions will be 
discussed. Finally, a summary and conclusions are drawn. 



II. EXPERIMENTAL METHOD AND DATA ANALYSIS 

The experiments were conducted at the Oslo Cyclotron 
Laboratory (OCL) using a 38-MeV ^He ion beam. Self- 
supporting targets of 206pb and ^''^Pb metallic foils, enriched 
to 99.8% and 99.9% and with thicknesses of 4.7 mg/cm^ and 
1 .4 mg/cm^, respectively, were used. The experiments ran for 
twelve days with a beam current of wO.6 nA. 
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The bombardment of ^He ions on the Pb targets opens a 
number of reaction channels, such as (^He, ^He'7), (^He, ay), 
(^He, xnay) and (^He, ^He'xny). The following reactions are 
analyzed in the present study: 

1. 208pb(3He,^He'7)^°^Pb 

2. 208pb(3He, a7)2''^Pb 

3. 206pb(3He, 3He'7)2''*^Pb 

4. 206pb(3He, ayf^^Vh 

Particle- 7 coincidences were measured for 205-208p^ using 
the CACTUS multidetector array ifTsIl . The charged particles 
were detected by eight coUimated AE— E type Si particle tele- 
scopes, placed at a distance of 5 cm from the target and mak- 
ing an angle of 45° with the beam line. The thicknesses of 
the AE and E detectors are «145 /xm and «1500 /xm, respec- 
tively. For the detection of 7-rays, 28 5"x5" Nal detectors 
were used, surrounding the particle telescopes and the target. 
The total efficiency of the Nal detectors is ^ 15% of An. 

The data analysis consists of four main steps: (i) making 
the particle- 7 coincidence matrix, (ii) unfolding the total 7- 
ray spectra, (iii) extracting the first generation 7-ray spectra 
and (iv) the factorization of the first generation 7-ray matrix 
into level density and 7-ray transmission coefficient. From the 
known Q values and reaction kinematics the ejectile energy 
can be transformed into initial excitation energy of the resid- 
ual nuclei. Using the particle-7 coincidence technique, each 
7-ray can be assigned to a cascade depopulating an excitation 
energy region in the residual nucleus. Thus, the particle-7 co- 
incidence measurements give a total 7-ray spectrum for each 
excitation energy bin. Each row of the coincidence matrix 
corresponds to a certain excitation energy (E) while each col- 
umn corresponds to a certain 7-ray energy (Ey). In 205,206p^ 
and ■^'^'»pb the excitation energy bins are chosen to be 240 
keV/channel and 220 keV/channel, respectively. 

The 7-ray spectra are corrected for the Nal response func- 
tion. The unfolding procedure of Ref. lfl6ll is employed for 
this purpose, which is based on the Compton-subtracting tech- 
nique that prevents additional spurious fluctuations in the un- 
folded spectrum. The reliability of the unfolding technique is 
tested for ^"^Pb 7-ray spectra, see Fig. [T] Here, the raw 7- 
ray spectrum from the excitation energy region E = 4.5 — 6.7 
MeV is compared with the spectrum that has been folded after 
unfolding. The good agreement with the raw 7-ray spectrum 
and the folded spectrum gives confidence in the employed un- 
folding technique. The set of unfolded 7-ray spectra are orga- 
nized in a {E,Ey) matrix, which comprises the energy distri- 
bution of all 7-rays from all decay cascades as a function of 
excitation energy. This matrix of total, unfolded 7-ray spectra 
is the basis for the next step of the Oslo method. 

The first generation (primary) 7-ray energy distribution is 
extracted by an iterative subtraction technique described in 
Ref. 1 17]. The basic assumption of the subtraction method 
is that the 7-decay pattern from any excitation energy bin is 
independent of whether the state is populated directly via scat- 
tering or neutron pick-up, or through 7-decay from a higher- 
lying excited state. 



The basic assumption may not be fulfilled if the direct re- 
action at lower excitation bins, whose 7-spectra are utilized to 
subtract higher-generation 7-rays from total specti'a at higher 
excitation energies, do not favor some levels within the ex- 
citation bin that are populated by 7-rays from above. Such 
situations may cause that some 7-rays are not fully subtracted 
from the total spectrum. 

The influence of a possible different selectivity of levels at 
one excitation energy in the direct reaction compared to 7- 
decay from higher lying levels, is expected to be most pro- 
nounced when only a few levels are present in the excitation 
bin. Thus, these considerations are increasingly important for 
nuclei in the vicinity of closed shells having low level den- 
sity. It is therefore necessary to be cautious when applying 
the Oslo method to the Pb nuclei, and to carefully check that 
the method gives reasonable results compared to other data. 

In Fig. |2]the application of the first generation 7-rays ex- 
traction technique is shown for ^"^Pb. The first generation 7- 
rays are extracted by subtracting the weighted sum of higher 
generation 7-rays from the total unfolded 7-ray spectrum. Fig- 
ure |2] shows the total unfolded 7-ray spectrum, the first gen- 
eration 7-ray spectrum, and the second and higher generation 
7-spectrum from E = 4.5 — 6.7 MeV for ^"^Pb. It is seen how 
the higher generation 7-rays are well separated from the total 
7-ray spectrum in the specified excitation energy region. 

in. NUCLEAR LEVEL DENSITIES 

A. Application of tlie Axel-Brink hypothesis 

The energy distribution of primary 7-rays emitted from a 
well-defined initial excitation energy provides information on 
the level density and the 7-strength function. First, the exper- 
imental primary 7-ray matrix P{E,Ey) is normalized for the 
various excitation energy bins E. This is done by summing 
P{E,Ey) over all 7-ray energies Ey, for each excitation energy 
bin E such that 

E 

PiE,Ey)^l. (1) 

77 cmin 

The entries of the first generation matrix are the probabilities 
P{E,Ey) that a 7-ray of energy Ey is emitted from excitation 
energy E. In accordance with Fermi's golden rule, the decay 
rate is proportional to the level density of final states and the 
square of the transition matrix element. Thus, entries of the 
primary 7-ray matrix P{E,Ey) in the statistical region are pro- 
portional to the level density p{E — Ey) and the 7-ray trans- 
mission coefficient 3^{Ey): 

P{E,Ey)c<p(E-Ey).'7{Ey). (2) 

In the factorization procedure of the two-dimensional primary 
7-ray matrix, a lower limit of the excitation energy region 
£■ 3.0 MeV was used for all the analyzed nuclei in order 
to exclude the main part of non-statistical transitions from 
our data. The normalized experimental primary 7-ray matrix 
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P{E,Ey) can be approximated by 



p{E-Ey)£r{Ey 



n^^ErPiE-Ey).Sr{Ey 



(3) 



The first trial function for p is assumed to be unity i.e., p = 1 
and the corresponding 5^ can be determined by Eq. ([3]). Then, 
a minimum is calculated for each data point of p and ^ si- 
multaneously. This procedure is repeated until a global least- 
square fit is achieved for all the data points of P{E,Ey). 
In Fig.|3]such a least fit is compared with the experimen- 

n S ^ 9 (Y7 

tal primary 7-ray matrix for the Pb( He, a/) Pb reaction 
at different initial excitation energies between E = 3.8 — 6.7 
MeV with energy bins of 220 keV. The calculated primary 
7-ray spectra (solid lines) are obtained by multiplying the ex- 
tracted p and as defined in Eq. (|3]i. 

The 7-ray transmission coefficient ^{Ey) in Eqs. (|2]l and 
(|3]l is independent of excitation energy according to the gen- 
eralized Axel-Brink hypothesis lfT8lll9ll . stating that collective 
excitations built on excited states have the same properties as 
those built on the ground state. The average temperature in 
the excitation energy region studied is below 1.0 MeV for the 
final levels. The temperature is believed to depend on the fi- 
nal excitation energy by T 0= ^JE^, which is a slowly varying 
function. Thus, we assume that using a constant temperature 
for the factorization of Eq. (|2]i is approximately valid in our 
excitation region. The error bars of the data points take into 
account only the statistical errors, see Ref. [5]. This means 
that neither possible shortcomings of the first-generation pro- 
cedure nor a weak dependence on the excitation energy in the 
transmission coefficient, which means that the Axel-Brink hy- 
pothesis is not fully valid, are included in the error bars. Keep- 
ing this in mind, the comparison shown in Fig. [3] indicates 
that the analyzing method works satisfactorily for the ^"^Pb 
nucleus. 

The multiplicative functions of Eq. give an infinite num- 
ber of solutions. It has been shown in Ref. |5] that if one so- 
lution is known for Eq. Q, then the product p • 5^ is invariant 
under the transformation 

p{E-Ey)=A exp [a {E -Ey)]p{E- Ey) , (4) 



.f{Ey) = Bexp{a Ey) .^{Ey). (5) 

Therefore, neither the slope nor the absolute value of p and 
can be determined directly from the iteration procedure. The 
free parameters A, B and a must be determined by experimen- 
tal data in order to normalize p and 3/^. 



B. Normalization of nuclear level density 

The parameters A and a in Eqs. (|4|i and Q are obtained 
by fitting the inferred data points to the known discrete lev- 
els 1I20I1 at low excitation energies and to the level density 
at the neutron separation energy S„. The level density at Sn 
can be deduced from the Fermi-gas expression ll2lll using the 
available proton or neutron-resonance spacing data ll22ll and 



assuming that positive and negative parities contribute equally 
to the level density at S„. For £ = capture (s-waves), the level 
density po becomes: 



Po(5„) 



2a2 



[(/, + l)exp(-(/, + l)V2(j2) 



2m-1 



(6) 



For £ = I capture (p-waves), the above equation becomes: 

PliSn) = — [(/, - l)exp(-(/, - 1)V2(t2) 

+ /,exp(-/2/2c7^) 

+ (/^ + l)exp(-(/, + l)V2cj2) 

+ (/,+2)exp(-(/,+2)V2cj2)]-i, (7) 

where Dq and Di are the average s- and p-wave resonance 
spacings. The parameter /, is the spin of the target nucleus. 
For spin /( = 0, the first two terms inside the bracket [• • •] of 
Eq. (|7]) should be omitted and for spin /, = 1 /2 and 1 only the 
first term should be omitted. The spin-cut off parameter a is 
defined in Ref. IH by 



(T^ = 0.0888A2/\/fl(£-£, 



(8) 



where A is the mass number, a is the level density parame- 
ter and Epair is the pairing correction parameter The pairing 
correction parameter is estimated following the description of 
Ref. 1123(1 . The spin distribution of levels at one excitation en- 
ergy is given by ll2lll : 



21+1 

2(72 



exp[-(/+l/2)2/2(j2], 



(9) 



which is normalized to 'E.ig{E,I) ~ 1. The spin assignments 
of 206,208p|^ jjjg excitation region E ^ 4-5 MeV are taken 
from Ref. [201. The average spin distribution of these data 
points are compared with the relative spin distribution deter- 
mined from Eq. (|9|l, as shown in Fig. ID Within the statistical 
uncertainty of the data points the agreement of the experimen- 
tal and theoretical spin distributions is very good and supports 
our adopted a in Eq. 

In Ref. 1 22] both s-wave and p-wave resonance spacings Dq 
and Di for the target 204,206, 207 nuclei are given at S„. The 
deduced level densities po and pi at S„ are listed in Table U 



Since the s-wave resonances in 



207.208 



Pb are more weakly 



populated than the p-wave resonances, the level densities de- 
termined using Dq are less reliable than those determined us- 
ing Di. Therefore, the pi level density has been used for nor- 
malizing our data. 

The parameters used in the analysis of 205-208p^ ^j.g ^^J^^_ 
marized in Table U The level density parameter a of the Pb 
nuclei is taken from Ref. (2^ using the Gilbert and Cameron 
approach ||2lll. Since the proton/neutron resonance spacing 
information for 206p^ jg unavailable, we estimate the level 
density at the neutron separation energy using the systemat- 
ics of Ref. ll2lll . Figure |3 demonstrates the estimation proce- 
dure, where the level densities of odd and even Pb nuclei are 



4 



shown. In ^"^^Pb and ^"^Pb the Fermi gas level densities at 
S„ are comparable with those deduced from the neutron res- 
onance data |22]. However, in ^*'^Pb and ^''^Pb the large dis- 
crepancy between the two level densities is obvious. In view 
of this discrepancy we have estimated an uncertainity of 80% 
for the adopted Fermi gas level density in ^"^Pb. 

The experimental level density is extracted up to ^ 2 MeV 
below S„ since the 7-rays below 2 MeV are omitted in the ex- 
traction procedure. In order to fill the gap between the data 
points and the level density at S„ deduced from resonance 
data, an interpolation is required. The Fermi gas level den- 
sity formula for all spins and parities 



exp{2\ 



12V2flV4f/5/4(j 



(10) 



Epair is the 



is employed for the interpolation. Here, U = E 
intrinsic excitation energy and Tj is a constant used to adjust 
Ppg to the experimental level density at S„ (values are given 
in Table IB. 

Figures |6] and [T] show the normalized level densities ex- 
tracted for 205,206p|3 ^jjj 207.208p|3^ respectively, up to 5 - 6 
MeV of excitation energy. The exponential increase of level 
density with excitation energy is evident from these figures, 
where the Oslo data points are shown as filled squares. The 
Fermi gas level densities (dashed lines) that are used for nor- 
malization at the highest energy points, are seen to describe 
the data points only to some extent. This feature is expected 
for nuclei near closed shells, having few nuclear levels. For 
the lead region, it is clear that the Fermi gas gives a poor de- 
scription below 4-5 MeV of excitation energies. 

By comparing our data with the discrete levels from spec- 
troscopic experiments |20], a very good resemblance is seen at 
lower excitation energies in all the nuclei. Local differences 
might be due to violation of the assumptions behind Eq. (|2|l 
for nuclei where level densities are low and large fluctuations 
of the 7-ray intensity are observed. In the doubly magic ^ospb 
the level density is low and the Oslo data agree nicely with the 
spectroscopic measurements [20] up to 5 MeV of excitation 
energy. The presence of a single unpaired neutron in ^"^Pb 
increases the level density compared to 208p^ ^j. ^ gjyen exci- 
tation energy. The Oslo data of ^''^Pb provide new information 
on level density above E ~A MeV. Similarly, in 205,206p^ 
information on the level densities are determined at higher ex- 
citation energies where spectroscopic methods fail to find lev- 
els. 

Figure |6] shows that the level densities of 205,206p^ 
smoother functions of excitation energy with an overall higher 
absolute value than of 207.208p^^ where characteristic struc- 
tures are prominent. This feature is expected due to the pres- 
ence of two and three neutron vacancies, which provide extra 
degrees of freedom for the nucleons to arrange themselves for 
a given excitation energy. 

The drop of the total level density for E <A MeV while go- 
ing from ^"^Pb to ^"^Pb, is interpreted as a shell closure effect. 
The extraction of level density in ^"^Pb does not provide any 
new information above the known levels. However, these re- 
sults give further confidence in the Oslo method and its appli- 
cability in this mass region, as our data show good agreement 



with the data found in literature. In summary, these results are 
gratifying and support the applicability of the Oslo method in 
closed-shell nuclei where the level densities are low. 



IV. THERMODYNAMICS 

The nuclear level density is closely related to the entropy 
S of the system at a given excitation energy E. In fact, the 
level density p{E) is directly proportional to the number of 
accessible levels at excitation energy E. In order to derive 
thermodynamic quantities for mesoscopic systems, it is 
common to use the micro-canonical or canonical ensemble. 
The micro-canonical ensemble is the most appropriate 
ensemble according to Ref. ||25|]. since the nucleus is con- 
sidered as an isolated system with well-defined energy. 
The micro-canonical ensemble theory will therefore be 
utilized for 205-208p^^ although temperature and heat ca- 
pacity may become negative due to fluctuations in the entropy. 

In the micro-canonical ensemble, the entropy S{E) is iden- 
tical to the partition function determined by the multiplicity 
of states £!(£') at excitation energy E which corresponds to 
the level density p{E). The entropy in the micro-canonical 
ensemble can be derived as 



S{E) = kB\nQ.{E) 

...in^ 
Po 

^kB\np{E)+So, 



(11) 



where Boltzmann's constant ks is set to unity for the sake of 
simplicity. The value of the normalization factor is ad- 
justed such that the third law of thermodynamics is fulfilled: 
5^0 for r ^ 0. Since the ground-state band of the even- 



even 20bp^ ^I'Sp^ nuclei has T = 0, the parameter used 
for 205-208p^ 0.34. From the entropy, one can derive the 
temperature by 



1 



dS 



T{E) dE' 



(12) 



In Figs. [8] and |9] the micro-canonical entropies of 205,206p^ 

and 207,208p,, 

are shown, respectively. The entropies of 
207,208p^ are seen to vary strongly with excitation energy. For 
^"^Pb, the first vibrational 3^ state appears at 2.6 MeV, be- 
fore the two quasi-particle regime which is seen to enter at 
£■ > 3.5 MeV. Figure [TO] shows that the difference in entropy 
between ^o^pb and ^o^Pb, AS = ^(^''^Pb) - ^(^ospb), is fluc- 
tuating between 0.3 — 2.0 in the excitation energy region of 
£■ = 3.0 — 5.7 MeV. This strong variation in entropy is due to 
the few available single-particle orbitals for the unpaired neu- 
trons below the closed shell. The strong fluctuation of entropy 
in 207,208p^ excitation energy makes it difficult to deter- 
mine other thermodynamic quantities such as temperature. 

In Fig. [To] the average difference in entropy between ^"^Pb 
and 206pb is shown to lie around 0.9 for the excitation en- 
ergy region between E = 2 — 5 MeV. The entropy difference 
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of ^'^^Pb and ^^°Ph is less fluctuating than that observed be- 
tween ^'^^Pb and ^"**Pb, thus indicating the departure of shell 
closure effects in 205,206p^ entropy difference between 
^''^Pb and ^"^Pb is about one half of the value observed in rare 
earth nuclei (AS ^ 2) |2dl. 

The micro-canonical temperature can be extracted from the 
entropy as given by Eq. (fT2l i. Since the level densities in 
207,208p|^ are low and show large fluctuations, the tempera- 
ture extraction for these nuclei is probably not reliable. The 
micro-canonical temperatures of ^"^Pb and ^''''Pb are shown in 
Figs.[TT]and [12] respectively. For ^"^Pb, the temperature ex- 
traction gives rise to large bump structures. These structures 
in the temperature arise from the differentiation of the entropy 
(see Eq. (fTSTi). and can be interpreted as the breaking of nu- 
cleon pairs. When nucleon pairs are broken, new degrees of 
freedom open up leading to an increase of p(£') and decrease 
in the temperature T{E). 

For ^''^Pb, the situation is not so clear since this nucleus has 
an unpaired neutron which gives a smoother entropy and thus 
less structures in the temperature. Bumps indicating the pair- 
breaking process might also be seen here, but one can also 
fit a constant temperature to the data points in Fig. [TT] By 
doing so, we get an average temperature of 0.9(1) MeV for 
0.4 < E < 5.0 MeV, which is in agreement with an average 
temperature of T = 0.81 (4) MeV 

In 205pb and ^'"^Pb one can observe the first drop in tem- 
perature at £■ w 2.0 MeV and E w 2.5 MeV, respectively. 
These excitation energies can be compared with the energy 
amount necessary to break a neutron Cooper pair, namely 2A„, 
where A„ is the neutron pair gap parameter The description of 
Ref. ImI gives 2A„ = 1.8 MeV for 205,206p^ observed lo- 
cations of pair-breaking lie at energies somewhat higher than 
2A„. This difference could be due to the large energy spac- 
ing between the single-particle orbitals. However, the second 
peak in ^o^pb and ^O^Pb is observed at £ w 3.8 and 3.7 MeV, 
respectively. These peaks may be due to proton-pair break- 
ing, neutron-pair breaking, or because the proton (or neutron) 
passes the energy gap between major shell gaps. 

V. GAMMA-RAY STRENGTH FUNCTIONS 

The 7-ray transmission coefficient in Eq. ^ is con- 
nected to the electromagnetic decay properties of the nucleus 
and is expressed as the sum of all the 7-ray strength functions 
fxL for transitions with electromagnetic character X and mul- 
tipolarity L, given by: 

^(£^)-27r£4L+i/xL(£r)> (13) 

XL 

where Ey is the transition energy, and the slope correction 
exp(a£'y) from Eq. (|4]i has already been included in In 
Fig. [13] £^ for ^"^Pb is shown, where the absolute normaliza- 
tion B remains to be determined. This normalization will be 
discussed in the following. 

As already mentioned in sect. Ill, the methodological dif- 
ficulties in the primary 7-ray extraction prevent the determi- 
nation of p for E >S„-2 MeV and ^ for Ey<2 MeV. The 



level densities were extrapolated with the Fermi gas level den- 
sity using Eq. ( fTOl i. However, the transmission coefficients are 
extrapolated using an exponential form, as shown in Fig.[T3] 

It is assumed that the main contributions to the function 
.^{Ey) are £1 and Ml 7-ray transitions in the statistical re- 
gion. Further, if one assumes that the number of accessible 
levels of positive and negative parity are equal for any energy 
and spin, one finds 

p{E-Ey,If,±Kf)=^-p{E-Ey,If). (14) 

Clearly, this assumption does not hold true for all excitation 
energy regions in the Pb nuclei considered here. However, by 
applying these assumptions, the experimental 7-ray transmis- 
sion coefficient of Eq. (fT3] ) yields 

B^{Ey) = 2K[fEx{Ey)+,fM\{Ey)]El, (15) 

where S is an unknown factor that gives the absolute normal- 
ization of the y-ray strength function. We might expect a po- 
tential error in the absolute normalization of the 7-ray strength 
functions in the Pb region due to equal-parity assumption used 
above. 

The average total radiative width {Yy) of levels with energy 
£, spin / and parity K is given by 12811 : 

xp{E-Ey,If,nf), (16) 

where the summation goes over spins of final levels If with 
parities Kf. Using the assumptions discussed above, one can 
determine the average total radiative width (Fy) for neutron 
resonances. For i = Q capture (s-waves), the populated spins 
are / = |/, ± 1/2|, where /, represents the spin of the target 
nucleus in the (n, 7) reaction. For £—1 capture (p-waves), the 
populated spins are/ = |/, ±1/2±1|. The parity is determined 
by the target parity tt, and ^ by tt = tt, (— 1 )^. The average total 
radiative width is described at £ = S„ by 

1 fS" 

{Ty(Sn,l)) =- / dEy 

■XB.^{Ey)p{S„-Ey) 

1 

X £ g{Sn-Ey,I + J). (17) 

J=-\ 

The factor B can thus be determined from the average total 
radiative width (Fy) of the compound states. Ref. |'22] pro- 
vides the experimental (Fy) for ^"^Pb, while the (Fy) in ^"^Pb 
is taken from Ref. 1 1 ] . 

In order to normalize the 7-ray strength function of ^ospb 
the above mentioned method is difficult to apply since the av- 
erage total radiative width at 5„ cannot be found in the litera- 
ture. However, there are discrete neutron resonance data |2§] 
for ^"^Pb on the E\ and Ml 7-ray strength functions at Ey = 
7.5 MeV. The Oslo data for ^ospb is thus normalized with the 
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data of Ref. 112911 . In addition, the (7,n) photonuclear cross- 
section data of Ref. fso'] has been used to confirm the normal- 
ization of the 7-ray strength function of ^"^^Pb. The following 
relation ll22ll can be used to deduce the strength function from 
the cross-section, asumming that dipole radiation is dominant: 



f{Ey 



a{Ey 



(18) 



Since ^"^Pb is unstable, the resonance spacings and (Ty) 
have not been measured for ^''^Pb. Therefore, we utilized the 
(7,n) cross-section data |30] to normalize the 7-ray strength 
function of ^"''Pb. The Oslo data are scaled to match the data 
of Ref. 1 30], giving the absolute normalization of the 7-ray 
strength function of ^"^Pb. The normalized 7-ray strength 

functions of 205,206p^ ^jjj 207,208p^ gj^^^jj pjg^ [g] 

and[T5] respectively. 

For ^"^Pb, the resonance data ||29ll for electric and magnetic 
dipole radiation at fiy = 7.1 MeV have been used to make a 
consistency check for the normalized 7-ray strength function. 
Moreover, the (7,n) data of Ref. Q for ^o^Pb at Ey > 5„ ai-e 
also drawn in Fig.[T5]for comparison. The agreement of these 
strength functions with our extracted data is gratifying, and 
indicates that the normalization procedure works well in this 
case, in spite of the questionable assumption of equal parity 
distribution described in Eq. (fl4] i. 

The (7,n) data for 206-208p^ Q display bumps and 

structures that increase in magnitude when approaching S„. 
Similar structures are present in our data below 5„. These in- 
termediate structures are observed in the Oslo data at 7-ray 
energies 5.6, 6.3 and 7.1 MeV in ^ospb. Similarly, in ^o^Pb 
these intermediate structures have been seen at 4.3, 5.0, 5.6 
and 6.3 MeV. For ^''^Pb and ^°^Pb these structures are ob- 
served at 7-ray energies of 4.3, 6.0 and 6.3 MeV. It is hard 
to interpret these structures without knowing the electromag- 
netic character and the multipolarity of transitions in this re- 
gion. However, it is likely that the bumps are due to enhanced 
single-particle transitions at certain 7-ray energies, reflecting 
the shell structure below 5„ . The Oslo method does not pro- 
vide any information on the multipolarity of these structures, 
so to investigate it further, it would be necessary to perform 
complementary experiments such as (n,77) measurements an- 
alyzed with the two-step cascade method [3 ill . 

The above extracted 7-ray strength functions depend on the 
normalization procedure chosen. The slope of the strength 
function is sensitive to the resonance data at S„, which have 
been taken from the literature. In 206pb the uncertainty of this 
data play a central role in the normalization of our data. The 
adopted values influence both the slope of the level density 
as well as the slope of the strength function. In ^"^Pb, the 
uncertainty in the value of a (see Eq. ( fT3] )) could be more 
than 80% due to this uncertainty. In addition, the absolute 
strength (B of Eq. (flST l) is uncertain by a factor of 2 - 3. In 
the case of ^"'^Pb, the oscillating shape of the strength function 
makes the absolute normalization between our data and (7,n) 
data f29ll30ll very difficult, as seen in Fig.fTSl 



VI. MODELS FOR El AND Ml TRANSITIONS 

A number of models describing the electric dipole 7-ray 
strength function have been described in |22]. The simplest 
of these is the Standard Lorentzian model (SLO), which de- 
scribes the £1 strength as a Lorentzian shape. The model is 
temperature independent and is given by: 



/sLo(£'r) ' 



1 



El 



{El 



{TElEy)^' 



(19) 



where the Lorentzian parameters (Jei, Eei and F^i are the 
peak cross section, energy and width of the giant dipole reso- 
nance, respectively, and are usually derived from photonuclear 
experiments. However, it has been shown |28, 29, 32, 33, 3^ 
that the SLO model overestimates the photonuclear data away 
from the giant dipole resonance (GDR) centroid in many nu- 
clei. 

One model that incorporates the temperature dependency 
of the 7-ray dipole strength is the Enhanced Generalized 
Lorentzian model (EGLO) [2Z1 . Considering spherical nuclei, 
EGLO is defined as 



/eGLO (Ey) 



1 



r<7Eir£ 



EyTi;{Ey,T) 



{El 



72 \2_ 



{EyYk{Eyj)Y 

0.7 ^ . ' L (20) 
Eex J 



where T is th e temperature of the final states determined by 
T — a/ {U I a). The energy and temperature dependent width 
TkiEy, T) is defined analytically as 



Tk{Ey,T)^K{Ey 



-{2nTf 



where 



K{Ey) = k+{1-k] 



-'7- 



-Eo 



Eei -Eq 



(21) 



(22) 



is an empirical function. Here we use £0 = 4.5 MeV and the 
enhancement factor K, given by ll35ll . 



1 



if A < 148, 



[1 +0.09(A-148)2exp(-0.18(A-148)) if A > 148, 

(23) 

for the Fermi gas model. These expressions are developed in 
the framework of the collisional damping model for Ey < Eei 
and hold for T < 2 MeV. 

The magnetic dipole Ml radiation is described by a 
Lorentzian based on the existence of a giant magnetic dipole 
resonance (GMDR) ||36ll . 



/mi (Ey) 



1 



GmiEvT^ 



Ml 



F2r2 
Ml 



(24) 
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where ffMi , Tmi , and Emi are the GMDR parameters deduced 
from the systematics given in I22ll . 

The contribution from isoscalar E2 transition strength is 
also included in the total 7-ray strength function and is de- 
scribed by Ii2ai 



fE2{Ey) = 



1 



r 



72 ^2 



"£2 



y +Ejr^2 



(25) 



The total model 7-ray strength functions, shown by solid lines 
in Figs. [T6l and [TTl are the sum of the EGLO El, Lorentzian 
Ml and E2 contributions, i.e.. 



/tot — /eglo + /mi 



-EyfEi- 



(26) 



The GEDR and GMDR parameters for 205-208p^ ^^^^^ 
from i22ll and are listed in Table II. The extracted normalized 
7-ray strength functions of 205-208p|^ plotted together with 
the models discussed above in Figs.ll6land[T7l The photonu- 
clear cross-section data of Ref. OOll have also been drawn for 
comparison. 

The intermediate structures at the tail of GEDR in the Oslo 
data are observed at different 7-ray energies in the analyzed Pb 
nuclei. These intermediate structures present in our data os- 
cillate between the above-mentioned models; however, none 
of them describe our data adequately for the whole energy 
region. In Figs. [16] and [IT] our extracted data points (filled 

squares) of 205,206p|, 207,208p|, (gjjj fQjiQ^ jjjg £QLO 

model for 7-ray energies > 4 MeV. For lower 7-energies the 
deviation between the theory and data points is obvious. This 
is not surprising, since one expects that the closed shell(s) will 
strongly influence the 7-decay, and thus preventing a good 
description of the 7-strength with smooth functions such as 
given by the EGLO model. 

The strength functions in 206-208p^ show an increase at 7- 
ray energies lower than ^3 MeV. This may be due to the possi- 
ble presence of strong non-statistical transitions, which are not 
correctly subtracted in the primary 7-ray extraction procedure 
and thus affecting the strength function at these 7-energies. 
Therefore, the enhanced 7-ray strength functions at low 7-ray 
energies are not conclusive. 



has also been investigated in the doubly magic 208p|3 nucleus 
and its neighboring 205-207p|^ nuclei. In contrast to the rare 
earth and mid-shell nuclei, these isotopes have low level den- 
sity so that one can expect strong non-statistical fluctuations 
of level density and 7-ray strength function. 

The level densities and 7-ray strength functions of 
205-208p|^ have been extracted simultaneously from the pri- 
mary 7-ray spectra. The comparison of our extracted level 
densities with spectroscopic measurements at low excitation 
energies gives good agreement within the experimental un- 
certainty. In 2'^'^Pb, the Oslo method could not give more in- 
formation on level density than previously known from dis- 
crete spectroscopy. However, the good agreement between 
our data and known levels in 208p^ indicates the robustness of 
the method for its use in closed shell nuclei. 

The level densities of 207,208p^ show significant step struc- 
tures, which is an interesting finding of this work. Such struc- 
tures are expected, partly due to strong shell effects at the 
Z = 82 and N ^ 126 shell closures, and partly due to the 
breaking of Cooper pairs. In 205,206p^^ these step structures 
are smoothed out as neutron valence holes come into play. 

From the extracted level densities the micro-canonical en- 
tropies of the respective nuclei are deduced. The average en- 
tropy difference AS between ^''^Pb and ^'^^Ph is found to be 
0.9; however, the entropy difference between ^'^^Pb and 208p^ 
varies violently with excitation energy and thus is difficult to 
use for finding other thermodynamic properties. The fluctua- 
tions in the entropy spectra are strongly enhanced in the tem- 
perature spectra. The average temperatures of the 205,206p|^ 
nuclei are found to be T w 1 .0 MeV. This is twice the temper- 
ature measured in rare earth nuclei. 

The 7-ray strength functions of '^"^^■^''^Pb show pronounced 
structures at energies below the neutron separation energy. 
However, the 7-ray strength function becomes smoother by 
the gradual opening of the neutron shell closure at = 126. 
Such structures have also been observed above the neutron 
threshold for (7,n) reactions in 206-208p^ nuclei. The multi- 
polarities of these intermediate resonances are unknown. The 
measured 7-ray strength functions of 205-208p|^ ^j.g poorly de- 
scribed by the SLO and EGLO models in the energy region 
considered here. This indicates that there is more interesting 
physics connected to the shell closure(s), which is revealed in 
the 7-ray strength functions of the Pb isotopes. 



VII. SUMMARY AND CONCLUSIONS 

The primary motivation of this work was to determine the 
level density and the 7-ray strength function in Pb isotopes 
near and at shell closure. The applicability of the Oslo method 
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Sn 


a Epair Do 


Po(5„) 






I? 




(MeV) (MeV^l) (MeV) (keV) (10^ MeV" 


') (keV) (lO-^MeV" 


') 


2"5pb 
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0.30 
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207pb 


6.738 
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TABLE II: GEDR and GMDR parameters, used for the EGLO and SLO model calculations. 



Nucleus 


Ee\ Tei CTei 


Emi 




(rr) 




(MeV) (MeV) (mb) 


(MeV) (MeV) (mb) 


(meV) 


2"5pbt 


13.59 3.85 514 


6.95 


4.0 1.16 


330 


206pb 


13.59 3.85 514 


6.94 


4.0 1.16 




207pb 


13.56 3.96 481 


6.93 


4.0 1.16 


455(50) 


208pb 


13.43 4.07 639 


6.92 


4.0 1.72 





t GEDR parameters taken from ^^Ph. 
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FIG. 1 : Unfolding of the /-ray spectrum of ^''^Pb measured for the excitation energy region £ = 4.5 — 6.7 MeV. Note the similarity of the raw 
and folded spectrum. 
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FIG. 2: Extraction of first generation y-rays (middle panel) from the total, unfolded y-ray spectrum (top panel) by subtracting the second and 
higher generation y-rays (bottom panel) in ^"^Pb between excitation energies of 4.5 and 6.7 MeV. 
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FIG. 3: Comparison of the normalized experimental primary y-ray spectra for the ^*'^Pb(^He,a)^*'^Pb reaction (data points) at various excitation 
energies and the fit (solid lines) using the factorization of Eq. l|2j. The excitation energy widths are 220 keV. 
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FIG. 4: Average experimental spin distribution (data points with error bars) of 206.208p^ compared to the spin distribution determined using 
Eq. ([9}. The error bars include only statistical errors. 
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FIG. 5: Estimation of tlie level density at neutron separation energy p(5„) for ^"*Pb. The filled triangles represent the deduced p{S„) for 
205,207-209p(j ^siJJg £ = I neutron resonance spacings [22.1 . The filled circles are the level densities determined by Eq. l llOt using parameters 
of Table U 
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FIG. 6: Normalization of experimental nuclear densities (filled squares) for Pb and Pb. The arrows indicate the fitting region of data 
points at low excitation energies with known levels (solid line). To bridge the gap between the level density at S„ (open square) and the data 
points, the Fermi gas level density (dashed Une) is used. The filled triangle and the open square symbols at S„ are based on £ = and 1 neutron 
resonance spacing data, respectively. 
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FIG. 7: Normalized nuclear densities (filled squares) of ^"^'^''**Pb (see caption of Fig.|6ll. 
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FIG. 8: Micro-canonical entropy for 205,206p|j 
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FIG. 9: Micro-canonical entropy for ^u/,2USp|j 
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FIG. 10: Difference in entropy between "^Pb and ""Pb (upper panel) and between Pb and Pb (lower panel). The entropy difference is 
given by AS = Seo — See, where Seo and 5ee is the entropy of the even-odd and even-even nucleus, respectively. 
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FIG. 11: Temperature as a function of excitation energy in Pb. 
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FIG. 12: Temperature as a function of excitation energy in Pb. 
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FIG. 14: Normalized y-ray strength functions for 2U5,206pjj ^ function of y-ray energy. The Oslo data (filled squares) of Pb are normalized 
to photonuclear data (open circles) taken from Ref. (30ll . 
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FIG. 15: Normalized y-ray strength functions for 207,208 -pj^g Oslo data (filled squares) of ^^^Pb is normalized using the data of Ref. Ot]- 
The ^"^Pb nucleus is normalized by scaling the Oslo data with the sum oi fE\+ fui (open square) taken from Ref. (29ll and the photonuclear 
data (open circles) of Ref. t3Q1 . 




FIG. 16: Normalized y-ray strength functions for Pb and " Pb. The Oslo data (filled squares) are compared with the SLO (dotted line) 
and EGLO (dashed line) model predictions. The Ml strength function is shown by the dashed-dotted line. The solid line is the total strength 
function which is the sum of the EGLO model, E2 and Ml. The photonuclear cross-section data from Ref. BQI are shown as open circles. 
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FIG. 17: Normalized y-ray strength functions for 207, 208 pj^ ^ function of y-ray energy (see caption of Fig.[76l). 



